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The state of a two-particle system is called 
entangled when its quantum mechanical wave 
function cannot be factorized in two single- 
particle wave functions. Entanglement leads 
to the strongest counter-intuitive feature of 
quantum mechanics, namely nonlocality. 1 ' 2 
Experimental realization of quantum entangle- 
ment is relatively easy for the case of photons; 
a pump photon can spontaneously split into 
a pair of entangled photons inside a nonlin- 
ear crystal. In this paper we combine quan- 
tum entanglement with nanostructured metal 
optics 3 in the form of optically thick metal 
films perforated with a periodic array of sub- 
wavelength holes. These arrays act as pho- 
tonic crystals 4-7 that may convert entangled 
photons into surface-plasmon waves, i.e., com- 
pressive charge density waves. We address the 
question whether the entanglement survives 
such a conversion. We find that, in princi- 
ple, optical excitation of the surface plasmon 
modes of a metal is a coherent process at the 
single-particle level. However, the quality of 
the plasmon-assisted entanglement is limited 
by spatial dispersion of the hole arrays. This 
spatial dispersion is due to the nonlocal dielec- 
tric response of a metal, which is particularly 
large in the plasmonic regime; it introduces 
"which way" labels, that may kill entangle- 
ment. 

The samples that we have studied are two 
1 mm x 1 mm "metal hole arrays" , each comprising 
a 200 nm thick gold layer perforated with a square 
grid of 200 nm diameter holes spaced with a 700 nm 
lattice constant; a typical SEM picture is shown in 
the insert of Fig. 1. They were produced at DIMES 
(see acknowledgements) by first defining, with elec- 



tron beam lithography, arrays of dielectric pillars on a 
0.5 mm thick glass substrate, subsequently evaporat- 
ing the gold layer onto a 2 nm thin titanium bonding 
layer, and finally etching away the pillars to leave the 
array of air holes. 

Transmission spectra of two hole arrays, measured 
with a spectrometer using normally incident white 
light, are shown in Fig. 1. One can clearly see the res- 
onances due to excitation of Surface Plasmons (SP) 
on either of the metal-dielectric boundaries. At these 
resonances the measured transmission can be orders 
of magnitude larger than the value obtained from clas- 
sical diffraction theory for subwavelength holes. 4,8 In 
a simple picture, the surprisingly large transmission 
is due to the coupling of photon to SP on one side of 
the metal, subsequent tunnelling of the SP through 
the holes to establish a SP at the other side, and final 
re-radiation into a photon. 5 Other prominent fea- 
tures in the spectra are the transmission minima as- 
sociated with Wood anomalies. 6 At this moment the 
theoretical description of the full transmission spec- 
trum is incomplete, but the role of the SP herein is 
well established. 5-7 The resonance used in our ex- 
periment extends from 800-830 nm and is associated 
with the (±1,±1) SP mode on the glass- metal inter- 
face. Peak transmissions of the two arrays are typ- 
ically 3% (dashed curve) and 5% (solid curve) at a 
wavelength of 813 nm; these values are much larger 
than the value of 0.55% obtained from classical the- 
ory. 8 The difference in transmission between the two 
nominally identical hole arrays is ascribed to produc- 
tion imperfections. 

As a mainly longitudinal, compressive electron 
density wave, a SP propagates along the polariza- 
tion direction of the incident light, following a certain 
dispersion relation. 9, 10 In order to confirm this for 
our samples we have measured polarization-resolved 




V\&velength [nm] 

Figure 1: Wavelength-dependent transmission of the 
two hole arrays used in the experiment. The dashed 
vertical line indicates the wavelength of 813 nm used in 
the entanglement experiment. Inset: Scanning Electron 
Microscope picture of part of a typical hole array. The 
scale bar is 2 /.im long. 

transmission spectra of one of the hole arrays for var- 
ious angles of incidence 8 inc of plane-wave radiation 
(see Fig. 2). Angle tuning is expected to shift the var- 
ious resonances in the transmission spectrum in dif- 
ferent ways. As our interest is in a SP that propagates 
along the (±1,±1) direction, we have varied the angle 
of incidence by tilting the samples along the diagonal 
(4> = 45°) axis of the square hole array. For a po- 
larization orthogonal to this tilting axis (Fig. 2a) the 
main peak at 810 nm splits and shifts for increasing 
6i nc ; for a polarization along the tilting axis (Fig. 2b) 
the main peak remains at the same spectral position. 
These results confirm the relation between polariza- 
tion and propagation and thus the 2D photonic band 
structure of the SPs. 

We generate polarization-entangled photons with 
the standard method of type II Spontaneous Paramet- 
ric Down- Conversion (SPDC) 11,12 depicted in Fig. 3. 
A 240 mW CW Krypton-ion laser beam at a wave- 
length of 406.7 nm is directed onto a 1 mm thick 
BBO nonlinear crystal, where the beam diameter is 
ss 0.50 mm (full width at 1/e 2 points). Inside the non- 
linear crystal a small fraction of the pump photons is 
down-converted into twin photons at the double wave- 
length; these are emitted along two intersecting cones. 
Polarization-entangled photon pairs are selected with 
pinholes at the crossings of these cones; the size of 
the pinholes (far-field diameter 5 mrad) was chosen as 
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Figure 2: Wavelength-dependent transmission of a hole 
array for various angles of incidence at the two polar- 
izations. The lowest curve is measured with the probe 
beam at normal incidence, while consecutive curves are 
measured at increasing angles of incidence (one degree 
steps) by subsequent tilting of the square hole array 
around a diagonal. These curves have been plotted with 
a subsequent vertical offset of 2%. The dashed lines in- 
dicate the 813 nm wavelength used in the entanglement 
experiments. The resonance at this wavelength shows 
a complicated splitting for a polarization orthogonal to 
the tilting axis (a), while it is stationary for polarization 
along this axis (b). 

a compromise between high yield and good quantum 
entanglement. Lenses L of 40 cm focal length produce 
a one-to-one intermediate image of the pumped area, 
which will be used (see below) to accommodate the 
hole arrays Al and A2. After passing through polariz- 
ers PI and P2 the entangled twin photons are focused 
through interference filters (10 nm bandwidth cen- 
tered at 813 nm) onto single-photon counting mod- 
ules (Perkin Elmer SPCM-AQR-14). The photode- 
tectors act as bucket detectors, i.e., they impose no 
further transverse mode selection (this is only done by 
the pinholes) . The signals from the two detectors are 
combined electronically in a coincidence circuit with 
a time window of 2 ns, to provide the two-photon 
coincidence rate. 

In a simplified picture, the generated polarization- 
entangled state is 

\9) = ±(\H 1 V i ) + e^xfla)) , (1) 
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Figure 3: Experimental setup, comprising a BBO gener- 
ating crystal, pinholes, 40 cm-focal length lenses L, po- 
larizers PI and P2, interference filters IF, single photon 
counters 1 and 2, and a coincidence circuit (2 ns). Beam 
walk-off is compensated for with the standard quantum 
eraser comprising a half-wave plate HWP at 45° and 
compensating crystals C with a thickness equal to half 
of that of the generating crystal. 11,12 The dotted ob- 
jects are present only in some experiments; they show 
the hole arrays Al and A2 at the image positions created 
by lenses L, or, alternatively, in the focus of the confocal 
telescope TEL (15 mm-focal length lenses). 

where the state represents the simultaneous 

emission and propagation of a iJ-polarized photon in 
beam 1 and a l/-polarized photon in beam 2. The 
H- and V^-directions are defined by the orthogonal 
birefringent axes of the generating BBO crystal and 
all spatial information is implicitly contained in the 
beam labelling. By tilting one of the compensating 
crystals (C in Fig. 3) the quantum phase 6 can be 
set. 

In a more accurate description, where the spatial 
information within the beams is properly accounted 
for, the two-photon state in a given transverse plane 
is 
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where the integration is over the transverse modes 
present in beam 1 and 2, labelled by wavevectors 
q\ (i=l,2). $hv (qi,q 2 ) is the two-photon prob- 
ability amplitude for finding a H-polarized pho- 
ton with wavevector q\ in combination with a V- 
polarizcd photon with wavevector q 2 , and vice versa 
for $vh (<fi> <7a)- These functions are similar to the 
one normally used, 13 ' 14 but now with polarization la- 
belling and taken specifically at frequency degeneracy. 
If the functions $hv and "IVh are identical, i.e., in the 



absence of polarization-dependent components, they 
can be factorized out and Eq. (2) reduces to the simple 
form of Eq. (1). If the functions are different, polar- 
ization and spatial information are intertwined and 
the degree of polarization entanglement is reduced. 
Mathematically, it is the overlap integral between the 
two ^-functions in Eq. (2) that determines the degree 
of entanglement, as it measures the spatial indistin- 
guishability of the polarization-entangled photons. 14 

In the absence of the hole arrays our setup pro- 
duces typically 3.2 x 10 4 coincidence counts per sec- 
ond, being w 25% of the single count rate. A measure 
for the purity of the quantum entangled state is the 
so-called visibility of the biphoton fringe. 11,12 This 
visibility was typically Vb° = 99.3% and V450 = 97.0% 
at polarization angles of 0° and 45°, respectively (see 
Table). The high value at 45° shows that the natu- 
ral preference of the generating BBO crystal for its 
birefringent axes (0° and 90°) was almost completely 
removed in our setup by the compensating crystals, 
which act as quantum erasers. 

Placement of the two hole arrays in the two "beam 
lines" of course leads to a dramatic reduction of single 
and coincidence counts. Coincidence count rates are 
typically 55 s _1 at the optimum setting of the detect- 
ing polarizers, which is consistent with the transmis- 
sions of the arrays given above. We again performed 
a measurement of the purity of the entangled state 
and found that the visibilities were now Vbo = 97.1% 
and V450 = 97.2%, respectively. In Fig. 4 the cor- 
responding fourth-order quantum interference fringes 
are shown for polarizer P2 fixed at 0° and 45°, respec- 
tively, and PI varying in steps of 10°. These measure- 
ments show that the quantum entanglement survives 
the transition from photon to SP and back. 

As a further demonstration of the SP-assisted 
quantum entanglement we performed a standard mea- 
surement of the so-called S-parameter, as described in 
ref., 11 on a singlet Bell state (6 — it in Eq. 1). This 
experiment, which took 16 runs of 100 seconds each, 
gave a value of \S\ — 2.71 ± 0.02, which is a vio- 
lation of the inequality of Clauser, Home, Shimony, 
and Holt, 15 by 35 standard deviations. 

Further experiments were performed on a setup 
with only a single array in one of the beams. The 
results for this case look very similar to those for the 
experiment with two arrays; again the entanglement 
was practically unaffected (see Table). This is to be 
expected as the two-photon wavefunction of Eq. (1) 
is perturbed by changes in either of its single-photon 
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Figure 4: Biphoton fringes, corresponding to fourth or- 
der quantum interference, as measured with the two hole 
arrays in place, for P2 fixed at 0° (solid curve) and 45° 
(dashed curve) and PI varying in steps of 10°. 

components; in principle a single array can kill all 
of the entanglement. The difference between the two 
single-array experiments, as reported in the table, is 
due to imperfections in array 2, which are also ob- 
servable in its (single-photon) polarization-dependent 
transmission. As the measurements using only ar- 
ray 2 gave results very similar to the situation with 
both arrays in place, these imperfections must be the 
key culprit for the somewhat limited visibilities in the 
two-array experiment. 

The most intriguing results of our single-array ex- 
periment are obtained when we focus the SPDC beam 
onto its hole array, using a confocal telescope (close 
to lens L) of two f = 15 mm lenses symmetrically po- 
sitioned around the array, as shown inside the dashed 
square in Fig. 3. Under these conditions we observe 
a drastic reduction of the degree of quantum entan- 
glement: when the intra-telescope focus has a nu- 
merical aperture of 0.13 rad we observe visibilities of 
V o = 73% and V Ab ° = 87% (see table). 

The observed reduction in visibility upon focusing 
can be naturally explained as a consequence of spa- 
tial dispersion in the hole arrays. Spatial dispersion 
is a little-known phenomenon that refers to the fact 
that the optical constants of a medium depend on 
the wavevector of the incident light as a consequence 
of the non-local relation between the optical excita- 
tion in a medium and the incident optical field. It 
can lead, surprisingly, to polarization anisotropy of 
a cubic crystal, since the cubic symmetry is broken 
due to the finite size of the optical wavevector 16 (As 
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Table 1: Biphoton fringe visibilities under various condi- 
tions. Here R stands for the measured coincidence count 
rate, Vo° and V450 for the measured visibility for one of 
the polarizers fixed at 0°, 45° respectively. 



an aside, we note that this effect has recently been 
demonstrated to have serious consequences for the 
use of high-precision crystalline optics for lithographic 
wafer production in the far ultraviolet 17 ). Spatial dis- 
persion is usually negligible as compared to frequency 
dispersion, i.e., the fact that optical constants de- 
pend on the frequency of the incident light. However, 
spatial dispersion can become significant in conduct- 
ing media, where the motion of free charge carriers 
causes non-locality over distances much greater than 
atomic dimensions. 16 SPs are an extreme case thereof 
as these entities are not at all local but instead propa- 
gate along the dielectric-metal interface at nearly the 
speed of light over distances of many optical wave- 
lengths. 3 ' 9 - 10 

As a result of this spatial dispersion, the near-field 
distribution of the photons that are reradiated at the 
backside of the hole array differs from the spatial pro- 
file of the "polarization-isotropic" incident photons. 
From symmetry arguments we expect a four-lobed 
profile that is oriented either along the array axes 
or along the diagonals, where the latter applies to 
the (±1,±1) SP that we probe, as sketched in Fig. 5. 
The unpolarized overlap region of these profiles cor- 
responds to the focused incident light, whereas the 
polarized elongations are a result of SP propagation. 
The presence of these polarized lobes introduces the 
possibility to distinguish the polarization of the pho- 
tons on the basis of their spatial near-field coordi- 
nates. This will then automatically kill part of the 
polarization entanglement, just as any "which way" 
information will do. A similar argument applies to 
the far-field, where the SP propagation creates a dif- 
ference between the two ^-functions in Eq. (2) after 
passage through the hole array. Note that the ob- 
served reduction in visibility is much stronger for Vqo 
than for V450 , contrary to what one generally finds 
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without using hole arrays. 11 ' 12 This observation is 
consistent with the fact that we excite SPs propagat- 
ing in the "diagonal" (45°) directions. 
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Figure 5: Schematic picture of the near-field at the 
back side of array Al when this is positioned inside the 
telescope. 

The reduction of the visibility due to this effect 
is expected to be significant when the propagation 
length I of the SP is comparable to the spot size of 
the focused incident beam, or rather its transverse 
coherence length. This is indeed the case for our ex- 
periment: based on the spectral width of the trans- 
mission peak at 810 nm, we estimate I w 4 /zm; the 
number of holes "covered" by the propagating SPs 
is thus about 10. In comparison, the spot size is 
w 8 ( = 0.5 mm x 15/800) and the transverse co- 
herence length (set by the 5 mrad far-field selection) is 
» 4 /zm. We note that the propagation length differs 
by an order of magnitude from the value I « 40 jum 
predicted for a smooth pure gold film; 3 we ascribe this 
difference to extra damping by the titanium bonding 
layer and the hole patterning. 

From a general perspective, the observed conser- 
vation of quantum entanglement for the conversion 
from photon — > surface plasmon — > photon is a first 
demonstration of the true quantum nature of SPs. 
All experiments on SPs so far do not probe their 
quantum nature, but only their semiclassical disper- 
sion relation c.q. wave nature. In this perspective, 
it is interesting to note that the true quantum na- 
ture of the photon was only established in 1977, in 
anti-bunching experiments. 18 ' 19 Furthermore, a sim- 
ple estimate shows that SPs are very macroscopic, in 
the sense that they involve some 10 10 electrons. Our 
experiment proves that this macroscopic nature does 
not impede the quantum behaviour of SPs, as they 



can be entangled with photons to yield the expected 
fourth-order interference. 

Based on the existing theory of spatial disper- 
sion 16 the observed reduction in visibility in the fo- 
cused setup is strongly linked to the square symme- 
try of the lattice and the (±1, ±1) character of the SP 
that we have used. In the near future we plan to study 
hole arrays based on other Bravais lattices and space 
groups, of which there are 5, respectively 17 different 
types in 2 dimensions. 20 A very interesting case will 
be the hexagonal lattice, where the spatial dispersion 
is expected to be independent of the orientation of 
the lattice. 

In conclusion, by addressing the topic of plasmon- 
assisted quantum entanglement we have combined 
two intriguing fields of research, being (i) quantum in- 
formation, and (ii) nanostructured metal optics. We 
hope that our work will stimulate studies of entangle- 
ment transfer to other condensed-matter degrees of 
freedom. 
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